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ABSTRACT: Using the path integral representation of the density matrix propagator 
of quantum Brownian motion, we derive its asymptotic form for times greater than the 

so-called localization time, (h/'jkT)? , where 7 is the dissipation and T the temperature of 
the thermal environment. The localization time is typically greater than the decoherence 
time, but much shorter than the relaxation time, We use this result to show that the 
reduced density operator rapidly evolves into a state which is approximately diagonal in a 
set of generalized coherent states. We thus reproduce, using a completely different method, 
a result we previously obtained using the quantum state diffusion picture (Piivs.Rev.D52, 
7294 (1995)). We also go beyond this earlier result, in that we derive an explicit expression 
for the weighting of each phase space localized state in the approximately diagonal density 
matrix, as a function of the initial state. For sufficiently long times it is equal to the 
Wigner function, and we confirm that the Wigner function is positive for times greater 
than the localization time (multiplied by a number of order 1). 



1. INTRODUCTION 



One of the simplest open systems that is amenable to straightforward analysis is the 
quantum Brownian motion model. This model consists of a non-relativistic point particle, 
possibly in a potential, coupled to a bath of harmonic oscillators in a thermal state. The 
quantum Brownian motion model has been used very extensively in studies of decoherence 
and emergent classicality (see for example, Refs. [1,2,3,4,5,6,7,8]). 

In the simplest case of a free particle of mass m in a high temperature bath, with 
negligible dissipation the master equation for the reduced density matrix p(x, y) of the 
point particle is, 



where a = Am^kT /h '. (More general forms of this equation, together with the derivations 
of it may be found in many places. See, for example, Refs. [2,9,5]). 

One of the most important properties of (1.1) (and also its more general forms) 
is that the density operator tends to become approximately diagonal in both position 
and momentum after a short time. This has been seen in numerical solutions and in 
the evolution of particular types of initial states for which analytic solution is possible 



A more precise demonstration of this statement was given in Ref.[16] by appealing to 
an alternative description of open systems known as the quantum state diffusion picture 
[17,18,19,20,21]. In that picture, the density operator p satisfying (1.1) is regarded as a 
mean over a distribution of pure state density operators, 



where M denotes the mean (defined below), with the pure states evolving according to a 
non-linear stochastic Langevin-Ito equation, which for the model of this paper is, 




(1.1) 



[10,6,11,7,12,8,13,14,15]. 




(1.2) 



|#) = ~H\^)dt ~\{L- (L)) 2 \m + (L — (L» |V) df(f) 



(1.3) 



2 



for the normalized state vector \ip), where H = p 2 /2m and L = ax. Here, d£ is a com- 
plex differential random variable representing a complex Wiener process. The linear and 
quadratic means are, 



The appeal of this picture is that the solutions to the stochastic equation (1.3) appear to 
describe the expected behaviour of an individual history of the system, and have been seen 
to correspond to single runs of laboratory experiments. For example, for the quantum 
Brownian motion model, the solutions tend to phase space localized states of constant 
width whose centres undergo classical Brownian motion [20,16,22,23,24]. The timescale of 
this process, the localization time, is at slowest of order (h/^kT)*, which is the timescale 
on which the thermal fluctuations overtake the quantum fluctuations [25,26,27]. For an 
initial superposition of localized states a distance I apart, localization initally proceeds on 
a much shorter timescale, of order h 2 /(£ 2 m'jkT) (which is often called the decoherence 
time [8,14]), thereafter going over to the slower timescale above. 

For us, the interesting feature of the quantum state diffusion picture is that it gives some 
useful information about the form of the density operator on time scales greater than the 
localization time. Given a set of localized phase space solutions \^pq), the density operator 
may be reconstructed via (1.2). This, it may be shown [16], may be written explicitly as 



Here, f(p, q, t) is a non-negative, normalized solution to the Fokker-Planck equation de- 
scribing the classical Brownian motion undergone by the centres of the stationary so- 
lutions. This is therefore an explicit, albeit indirect, demonstration of the approach to 
approximately phase space diagonal form on short time scales. 

The above demonstration was described by us in detail in Ref.[16]. However, we were 
not able to deduce an explicit form for the function f(p, g, t) using the quantum state 
diffuion picture. That is, we know that it is a solution to the Fokker-Planck equation, 
but it was not clear how to pick out the particular solution corresponding to a particular 



M[didC] = dt, M[«] = 0, M[d£] = 



(1.4) 




(1.5) 
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initial density operator. Intuitively, it is clear that f(p, q, t) is something like the Wigner 
function of the initial state, coarse-grained sufficiently to make it positive, evolved forwards 
in time, and with the interference terms thrown away. We would like to be able to show 
this explicitly. 

The aim of the present paper is to derive the form (1.5) for times greater than the 
localization time directly from the path integral representation of the density matrix prop- 
agator corresponding to (1.1), without using the quantum state diffusion picture. As we 
shall see, this derivation has the advantage that it gives an explicit expression for f(p, q, t). 
In particular, we shall show that f(p,q,t) coincides with the Wigner function Wt(p,q) of 
the density operator at time t, for sufficiently large times. 



2. THE DENSITY MATRIX PROPAGATOR 



The solution to the master equation (1.1) may be written in terms of the propagator, 



J, 



Pt{x,y) = / dx dy J(x,y,t\x ,y ,0) po{x ,y ) 



(2.1) 



(see, for example, Refs.[2,27] for further details of the quantum Brownian motion model). 
The propagator may be given in general by a path integral expression, which for the 
particular case considered here is 

J(x f ,y f ,t\x ,yo,0) = J VxVy exp J dt(x 2 - y 2 ) - y J dt(x - yf^j (2.2) 

This is readily evaluated, with the result, 



ira 



J (xf, y f , t\xo, yo, 0) = exp ^— y(xj - x ) - {yj - y ) 

(xf-yj) 2 + {x f - y f)(xQ - yo) + (xq - y Q f ) (2.3) 



o?t 

(For convenience we will ignore prefactors in what follows. They may be recovered where 
required by appropriate normalizations.) 
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The main result of the present paper comes from the simple observation that the real 
part of the exponent in the path integral (2.2) may be written 

^(-*J* (x - v ?) = Jv*^f« (x - t? -*f (2.4) 

The path integral representation of the propagator may therefore be written, 

J(xf,yf,t\xo,y ,0) = J Vx K x (xf, t\xQ, 0) K^(yf,t\y ,0) (2.5) 

where 

Kx(xf, t\xo, 0) = J Vx exp J dt x 2 - a 2 J dt (x - xf^j (2.6) 

For a pure initial state, po(x, y) = ^Q(x)^Q(y), the density operator at time t may therefore 
be written, 

pt(x,y) = Jvx *x(x,t)*Uy,t) ( 2 - 7 ) 
where the (unnormalized) wave function ty x is given by 

*x(x f ,t) = J dx K x (x f ,t\x ,0) ^ (x ) (2.8) 

(Wave functions of this type often appear in discussions of systems undergoing continuous 
measurement [28,29,30,31].) 

Our strategy is to first evaluate the quantity K x , examine is asymptotic form for times 
greater than the localization time, and then use it to reconstruct the density matrix prop- 
agator, J. The reason we expect this to yield the desired result is that up to normalization 
factors and ignoring the fact that x is real not complex, Eq.(2.8) is essentially the solu- 
tion to the Langevin-Ito equation, (1.3), so the phase space localization effect should be 
visible in its long time limit. Moreover, Eq.(2.7) is the analogue of (1.2) or (1.5), so by 
reorganizing the functional integral over x(t), we might reasonably expect to derive (1.5). 

The path integral (2.6) is essentially the same as that for a harmonic oscillator coupled 
to an external source, with the complication that the frequency is complex. The path 
integral is therefore readily carried out (see Ref.[32], for example), with the result, 

(i o o i 
-Cl{Xf+X Q ) + -C 2 XfXQ 

+c$Xf + c^xq + C5 — a 2 J ds x 2 (s)^j (2.9) 



where, 



moo cos u>t 

C\ = : 2.10 

2smoot v ' 

(2.11) 

smoot 
2a 2 rt 

C3 



smoot Jo 

2 rt 



/ ds x(s) sinoos (2-12) 
JO 

2a L [ l 

ca = / ds x(s) smuj(t-s) (2.13) 

smut Jo 

ck = ds ds' x(s)x(s') sincuft — s) sinks' (2.14) 

moosinoot Jo Jo 



Here oo = a(l — i) and 



%a 2 \ 2 f^kT 



I . i 
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The timescale of evolution according to (2.9) is therefore a -1 , which coincides with 
the localization time discussed in Ref.[16]. The asymptotic properties of K x are now easily 
seen. As t^oo, C2^0 and ci^^ma(l+i) like e~ at . Since C2^0, the propagator K x factors 
into a product of functions of xq and xj. The wave function (2.8) therefore "forgets" its 
initial conditions and becomes proportional to a Gaussian of the form 

exp ^c\xj + c^xj^j (2.16) 

on a timescale a . This is in complete agreement with the quantum state diffusion picture 
analysis of Refs. [20,16]. 

Now introduce 

q = Re C3, p = % (Re C3 + Im C3) (2-17) 
ma 

Then the Gaussian may be written 

(i o \ A ... , 2 z_ ma^ .,0 

expl-cix +c 3 x J =expl - — (1 - i)(x - q) + -px + —{I - i)q 

= (x\*pq) eWQ-*)? (2.18) 

The propagator K x therefore has the form 

K x (x f ,t\x ,0) =N (x f \*pq) e^ 1 "*)* 

x exp ^— ciXq + C4X0 + C5 — a 2 y <is £ 2 (s)^ (2-19) 



The generalized coherent states l^pq) depend on x(t) only through p and q, which are 
functionals of x(t). They are close to minimal uncertainty states, satisfying, ApAq = h/V2 
[20,16]. 

The desired form of the propagator is now obtained by inserting (2.19) in (2.5), but 
reorganizing the functional integral over x(t) into ordinary integrations over p and q and 
functional integrations over remaining parts of x(t). This may be achieved by writing the 
functional integral over x(t) as 

J Vx = J dpdq J Vx S(p -p)S(q-q) (2.20) 

with p and q given in terms of x by (2.17). We thus obtain, 

// met 2 

dpdq / Vx S(p-p) S(q-q) (x f \^ pq ) {^ pq \yj) e~ q 

x exp (^ciXq - jfXvl + c a xq + clyoj 

x exp ^c 5 + eg - 2a 2 J ds x 2 (s)^j (2.21) 

This may be written, 

J{xf,yf,t\x ,yo,0) = J dpdq f(p, q,t\xQ, yo) (xf\^pq) {®pq\yf) ( 2 - 22 ) 



where 



f(p, q,t\xQ, y ) = J Vx 8(p - p) 5(q - q) e"^^ 



x exp ( ^ciXq - jfiyl + c a xq + c|y 

<-t 



x exp ^c 5 + c% - 2a 2 J ds x 2 (s)^j (2.23) 

We have clearly cast the result in the desired form. Folding an arbitrary initial state into 
the expression for the density matrix propagator (2.22), we obtain an expression of the 
desired form (1.5), where f(p, q, t) is given explicitly by, 

fiP, q,t) = J dx dy f(p, q, t\xo, yo) po(x , y ) (2.24) 

This is our first main result. 



3. THE PHASE SPACE DISTRIBUTION FUNCTION 



It remains to evaluate the path integral expression (2.23). To do this first notice that 
(2.23) may be written 

f(p, q, t\x Q , yo) = exp ^cixg - % -c\y\ + ^Q^j j dkdk' e i kp+ i k ' q 
x J Vx exp {^j^V ~ ^k'q + c^xq + c|?/o^ 

x exp ^c 5 + cl - 2a 2 J ds £ 2 (s)^ (3.1) 

The functional integral over x is a Gaussian, since C5 is quadratic in x and p, q and C4 are 
linear in x, but it involves inverting the functional matrix contain in the last exponential 
in (3.1), which does not look particularly easy. However, we are saved from having to do 
this calculation by the following observation. From Eq.(2.5) and Eq.(2.9) (for at » 1), 
we see that 

J(xf,yf,t\x ,y ,0) =exp (^c^x 2 + Xq) - ^c\(y 2 + 

x J Vx exp (c 3 x f + c%yf + c 4 x + c\yo) 

x exp ^c 5 + eg - 2a 2 J ds x 2 (s)^j (3.2) 

This functional integral over x in this expression is very similar in form to (3.1) but we 
already know what the answer is: it is Eq.(2.3). In particular, equating (3.2) and (2.3), 
we obtain 

J Vx exp (c^x f + c^yj +C4XQ + c\yo + C5 + C5 - 2a 2 J ds x 2 {s) 

(xf -xq) 2 - {yj -yo) 2 ] 



exp 



/ im 
a 2 t 



x 



(xf - y f ) 2 + (x f - yf)(x - y ) + (x - y ) 2 ^ 
exp ("^i(x 2 + x 2 Q ) + l -c\{y) + y 2 ^j (3.3) 



Now the point is that the formula (3.3) is true for arbitrary Xf,yf. In particular, using 
(2.17), we see that 

c?,Xf + c* 3 yf = — [ X f + yf- i(x f - yf)] q + -p(x f - y f ) (3.4) 
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Hence the functional integral (3.3) is exactly the same as the one appearing in (3.1) if, in 
(3.3), we make the substitutions 

(x f - y f ) -> - k, ™ [x f + y f - i(x f - y f )] -> - % -k' (3.5) 

Inverting for xj and yf, we therefore find that the functional integral over x(t) in (3.1) is 
equal to the right-hand side of (3.3) with 



Using this result, and changing variables from k' to K = k + k' /ma in (3.1), we obtain 

,/ \ f ma 2 . .m^O^O a<2t c 2\ 
f[P,q,t\x ,yo)=expl—q +t — — — £ J 

f n <r/ ( ( a2 t ma\ o ma r/ 2 (ma m \ , \ 
x / dkdK exp - — )k - — K z + — — ) fcif 

( % ( mXo ha 2 t \\ 
X 6XP V h V ~ maq H 1 Z_ 6 _ j j 

x exp ^mag + ( 3 - 8 ) 

where = ^(^0 + 2/0 ); £o = x ~~ V0- This may now be evaluated. 

An alternative way of writing (3.8) is to carry out the same steps, but to change 
variables in (3.1) from k, k to Xf, yf, with the formal result, 

f(p,q,t\x ,yo) = J dxfdy f exp {^i 1 ~ i)( x f ~ Q) 2 ~ ^P x .f^j 

-2fi( 1 + i )(yf-Q) + ftPyf) J(xf,yf,t\x , yo,Q) (3.9) 
Folding in the initial state via (2.24), we obtain, 

f ( ma \2 i ^ 

f{P,q,t)= dxfdyfex.pl — (l-i)(xf-q) - -px / I 

/ ma . ... o * \ / \ m r.\ 

x exp I — (1 + - + jVVf j Pt{xf,yf) (3.10) 

which has the appearance of a formal inversion of the relation (1.5). 

Because the coordinate transformation (3.6), (3.7) is complex some attention to the 
integration contour is necessary. In particular, k and k' are integrated along the real axis, 
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therefore Xf + yj is integrated along a purely imaginary contour and xj — yf along a real 
contour. More precisely, let X = \{xj + yj) and £ = xj — yj. Then (3.9) becomes 

/(p, ? , * i*o, »o) = £2 rfx /T ^ exp (t ( (x _ q)2 + ?) ~ ^ ma ^ x - ? ) - 

x J(X + |,X-|,t|x , 2/0,0) (3.11) 

Explicitly, this integral reads, 

/ioc r+oo f ma ( i i \ 

dX J d£ exp f — ((I -q) + — J - -ma£(X - q) - j^pU 

x exp {^{X - X )(C ~ Co) " + ££0 + £ 2 )) (3-12) 

where and £o defined in the same way as X and £. The X integral will clearly converge 
since the contour is along the imaginary axis, and the £ integral will converge for sufficiently 
large at. 

Letting X^X + q, the integral over X is readily carried out, with the result 

f(p,q,t\x ,yo)= J #exp (-Lp£ + ^(q - X )(£ - &) " + «0 + Co')) 

(3.13) 



x exp 



ma 

In" 



e 2 +fe- (e " eo)x2 



at 



The integral over £ may now be evaluated but it is not necessary to do this, since the form 
of the answer is now clear. For at » 1, the terms in the second exponential are negligible 
compared to the similiar terms in the first. Furthermore, the remaining terms have the 
form of the Wigner transform of the propagator [27,33]. We thus have the simple result, 

f(p, q, t\x , yo) » j ' di e'irt jfa ^ t \ Xo , £ , 0) (3.14) 

Attaching an arbitrary initial density matrix, it then follows from (2.24) that 

f(p,q,t) « J dte-iPt pt ( q+ l -^q- l -t) 

= W t (p,q) (3.15) 

That is, for at » 1, f(p,q, t) is the Wigner function of the density operator at time t. 
This is the second main result of this paper. 
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From any of the above representations of f(p, q, t) (other than (3.14)), or from Ref.[16], 
it is straighforward to show that f(p, q, t) obeys the Fokker-Planck equation, 

df _ _P df kT d 2 f + i d 2 f + h d 2 f 

dt m dq dp 2 dpdq 2m dq 2 

As we have seen, f(p,q,t) approaches the Wigner function Wt(p, q) for at » 1, which 
obeys the Fokker-Planck equation of classical Brownian motion: 

dW pdW 1rT1 d 2 W , n . 

— = --^— + 2m>yKr— T 3.17 
at m dq op 1 

What happens is that the last two terms in Eq.(3.16) become negligible for large at, as 

may be seen by studying the Wigner function propagator (below). 



4. THE POSITIVITY OF THE WIGNER FUNCTION 

We have shown that the density operator approaches the form (1.5), where f(p,q,t) 
is given by the Wigner function. However, f(p,q,t) is by construction positive, yet the 
Wigner function is not guaranteed to be positive in general [33]. What happens is that 
the Wigner function becomes strictly non-negative after a period of time, under evolution 
according to (the Wigner transform of) Eq.(l.l), as we now show. 

The Wigner transform of the relation (2.1) yields, 

W t (p, q) = J dpodqo K(p, q, t\p , q , 0) W (p , q ) (4.1) 

where K(p, q, t\pQ, qq? 0) is the Wigner function propagator, and is given by [27], 

K{p, q, t\p , q , 0) = exp ^-^(p - po) 2 - v {^q - qo - + a(p - p ) (q - q - 

(4.2) 

where, introducing D = 2rwykT, 

1 3m 3m 

^=D? "=dF> a= D? (43) 

It is well-known that the Wigner function may take negative values only through 
oscillations in 7i-sized regions of phase space, and that it may be rendered positive by 
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coarse-graining of such a region. Considered for example, the smeared Wigner function 

W H (P, q) = 2j dp'dq' exp f ^J^— ~ ^f-) W P^ ' ^ ^ 

This object is called the Husimi function [34]. It is equal to the expectation value of the 
corresponding density operator in a coherent state (of position width cr q ), (p, q\p\p, q), so 
is non-negative. 

Loosely speaking, what happens during time evolution according to (4.1), is that, after 
a certain amount of time, the propagator effectively smears the Wigner function over a 
region of phase space greater than h, and it becomes positive, in the manner of (4.4). We 
will now show this explicitly. 

Letting po~^po + p and qo^qo + q — Pot/m in (4.1) yields, 

W t (p, q) = J dp dqo exp [~PPq - vq\ + ap q ^ W (p +p,qo + q- ^) (4.5) 

The further transformation po^po + j^Qo yields, 

W t (p, q) = J d PO d qo exp (-PPq - Pq^j W (po + -^q +p,qo + q- ^) (4.6) 

where (3 = — . These two transformations are canonical, and therefore the trans- 
formed Wigner function appearing in the integrand of (4.6) is still the Wigner function of 
some state (unitarily related to the original one). Hence, 

W t (p, q) = J d PO d qo exp (-PPo - W pq (po, q ) (4.7) 

for some Wigner function W pq depending on p, q. This may now be recast as the smearing 
of a Husimi function: 

J dp dq exp (- (P '~ 2 ^° - dlo) Ww(P0, ») (4.8) 
The integral over po, qo is a Husimi function with a 2 = 1/(2(3). Hence Wt{p,q) > 
provided the integral over p' in (4.8) exists. This will be the case if p^ 1 > h 2 (3, that is, if 

Va\V ft \5 



4 >IVJ (W) (4 ' 9) 

The Wigner function will therefore be non-negative for times greater than the localization 
time (mulitplied by a number of order 1). 
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5. DISCUSSION 



We have shown that for times greater than the localization time, (h/'jkT) 2 , the density 
operator satisfying (1.1) approaches the form 



where Wt(p, q) is the Wigner function and the \^ pq ) are close to minimum uncertainty 
generalized coherent states. The Wigner function is strictly non-negative for times greater 
than the localization time (times a number of order 1). 

Diosi has also discussed the possibility of the phase space diagonal form (1.5) under 
evolution according to the master equation (1.1) [35]. His method was very different to 
ours, in that he used the properties of the coherent states to regard (1.5) as an expansion 
of the density operator. He found that such an expansion is possible for times greater than 
the localization time, times a number of order 1, in tune with our results. 

An advantage of deriving (5.1) using path integral methods, rather than quantum state 
diffusion, is that it yields and explicit expression for the phase space distribution function 
f(p,q,t). Another advantage is that it is not obviously restricted to Markovian master 
equations. The quantum state diffusion picture, in its current state of development, exists 
only for systems described by a Markovian master equation. It may exist in the non- 
Markovian case, but is yet to be developed. The exact propagator for quantum Brownian 
motion, for quadratic potentials, can be given in terms of a path integral [5], and is (mildly) 
non-Markovian. Since the method described here utitlizes path integrals, rather than the 
quantum state diffusion picture, there is a chance that our method may be valid in the 
non-Markovian case also, but this is still to be investigated. 

We have concentrated in this paper on the simplest possible model of quantum Brown- 
ian motion: the free particle in a high temperature environment with negligible dissipation 
7. It is clear, however, that remaining in the context of a Markovian master equation, it 
would be straighforward (although perhaps tedious) to extend our considerations to the 
case of a harmonic oscillator with non-trivial dissipation. In the quantum state diffusion 




(5.1) 
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picture analysis, this case was covered in Ref.[16] and we expect the path integral treatment 
of the present paper to yield comparable results. 

It is perhaps enlightening to comment on the various timescales involved in a more 
general quantum Brownian motion model, and sketch the expected general physical picture 
part of which is described by the results of this paper. 

In this paper, we have largely been concerned with the localization time, {Uj^kT)^ 1 
which is the timescale on which an arbitrary initial density operator approaches the form 
(1.5). The nomenclature "localization time" comes from the quantum state diffusion pic- 
ture, which was the picture first used to derive some of the results described in this paper. 
It is so named because it is the time scale on which an arbitrary initial wave function 
becomes localized in phase space under evolution according to Eq.(1.3) [16,20]. 

Also relevant is the decoherence time, h 2 /(£ 2 m-fkT), which is the timescale on which 
the off-diagonal terms of the density matrix are suppressed (in the position representation) 
[14]. The decoherence time necessarily involves a length scale i, which comes from the 
initial state. It could, for example, be the separation of a superposition of localized wave 
packets, and the decoherence time is then the time scale on which the interference between 
these packets is suppressed. 

If one is interested in emergent classicality for macroscopic systems, it is appropriate 
to choose values order 1 in c.g.s. units, for i, T, m and 7. The decoherence time is then 
typically much shorter than the localization time. This is in turn typically much shorter 
than the relaxation time, 7 _1 , which is the time scale on which the system approaches 
thermal equilibrium (when this is possible). 

Hence the general picture we have is as follows. Suppose the initial state of the system 
is a superposition of localized wave packets. Then the interference terms between these 
wave packets is destroyed on the decoherence timescale. After a few localization times 
the density matrix subsequently approaches the phase space diagonal form (1.5). After a 
much longer time of order the relaxation time, the system reaches thermal equilibrium. 
Discussions of emergent classicality usually concern times between the decoherence time 
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and the relaxation time, and it is this range of time which has been the primary concern 
of this paper. 
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